Measurement of the Decay Kl — > n°e + e 7 
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We report on a new measurement of the branching ratio B(Kl — > Tv°e + e~'y) using the KTeV 
detector. This analysis uses the full KTeV data set collected from 1997 to 2000. We reconstruct 139 
events over a background of 14, which results in B(Kl — > ir°e + e~j) — (1.62 ±0.14 s tat ± 0.09 sy st) x 
10 -8 . This result supersedes the earlier KTeV measurement of this branching ratio. 



PACS numbers: 13.20.Eb, 11.30.Er, 12.39.Fe, 13.40.Gp 
I. INTRODUCTION 



The decay Kl — > 7r°e + e~7 can be used to study the 
low-energy dynamics of neutral K mesons. In partic- 
ular, this decay is an important check of Chiral Per- 
turbation Theory (ChPT), which has been used to de- 
scribe kaon decays in which long distance effects domi- 
nate. Up to 0{p A ) in chiral perturbation theory, there 
are no free parameters and one predicts the branching 
ratio to be approximately 1.0 x 1CK 8 1]. In the related 
decay K L — > 7r°77, the 0(p 4 ) calculation was found to 
underestimate the measured branching ratio by a factor 
of three[l, i, i, i, @. To match the K L -> tt°77 data 
it was found necessary to extend the calculation to in- 
clude 0(p 6 ) terms while introducing a free parameter, 
ay. This parameter characterizes the contributions from 
vector meson exchange terms. 0- The addition of these 
effects to the Kl — > 7r°e + e _ 7 calculation results in an 
increase in the branching ratio to 2.4 x 1CP 8 , approx- 
imately twice the 0(p A ) calculation. Our new branch- 
ing ratio measurement can distinguish between the 0(p i ) 
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and 0(p 6 ) predictions. 

Two previous experimental results have been reported 
on this decay mode[l, |jj |T(3]. The most recent mea- 
surement comes from the KTeV 1997 data set and is 
based on 48 events with a background of 3.6 ± 1.1 
events. That measurement yielded B(Kl — * 7r°e + e~7) 
= (2.34 ±0.35 ±0.13) x 10~ 8 , where the rate was normal- 
ized to a now-obsolete value of B(Kl — > 7r°7r°). Using the 
latest measurement of B(K L ->■ 7r°7r°)[ll the 1997 
measurement can be rescaled to (2.17±0.32±0.12) x 10 -8 . 
We report here on a new measurement of this mode from 
the KTeV experiment. For this analysis, we use improved 
techniques to reanalyze the 1997 data set and combine it 
with a new measurement from the 1999 data set. 

The Kl — > 7r°e + e~7 decay can also be used to help 
understand the CP violating decay, Kl — > 7r°e + e~. The 
Kl — > 7r°e + e _ decay contains both CP violating and CP 
conserving amplitudes. Since the Kl — ► 7r°e + e _ 7 decay 
proceeds through a two photon intermediate state, it can 
be used to determine the CP conserving components in 
Kl — > 7r°e + e~, and thus allow one to determine the CP 
violating contribution in Kl — > 7r°e + e~. Also, because 
the rate for Kl — > 7r°e + e~7 is orders of magnitude higher 
than the rate for Kl — ► 7r°e + e _ , a better understanding 
of the Kl — > 7r°e + e _ 7 decay will help to reduce the 
backgrounds to the Kl — * 7r°e + e~ decay. 
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FIG. 1: Schematic of the KTeV detector. 



II. THE KTEV DETECTOR 

We collect Kl — > 7r°e + e~7 events using the KTeV 
detector located at Fermilab. The data analyzed were 
taken during the 1997 and 1999 rare decay running pe- 
riods and comprised 2.9 x 10 11 and 3.8 xlO 11 kaon de- 
cays, respectively. The KTeV experiment employed two 
different configurations during its operation. The E799 
configuration was used for this measurement and was op- 
timized for reconstructing rare kaon decays. 

In the KTeV experiment [l3| neutral kaons are pro- 
duced in interactions of 800 GeV/c protons with a beryl- 
lium oxide target. The resulting particles pass through a 
series of collimators to produce two nearly parallel beams. 
The beams also pass through lead and beryllium ab- 
sorbers to reduce the fraction of photons and neutrons 
in each beam. Charged particles are removed from the 
beams by sweeping magnets located downstream of the 
collimators. The vacuum decay volume begins approxi- 
mately 94 meters downstream of the target, far enough 
so that the majority of the K$ mesons have decayed, 
and extends to approximately 159 meters from the tar- 
get. The decay volume is surrounded by photon veto 
detectors that reject photons at angles greater than 100 
milliradians. 

The most critical detector elements for this analysis 
are a charged particle spectrometer and a pure Csl elec- 
tromagnetic calorimeter [14j. The KTeV spectrometer is 
used for reconstructing charged tracks. This spectrome- 
ter consists of four planes of drift chambers; two located 
upstream and two downstream of an analyzing magnet 
with a transverse momentum kick of 0.205 GeV/c. Each 
drift chamber contains four planes of wires; two to mea- 
sure the horizontal position and two to measure the ver- 
tical position with a precision of approximately 100/im. 
During data taking in 1999 the momentum kick was 
reduced to 0.150 GeV/c to increase the acceptance for 
multi-track events. 

The Csl calorimeter is composed of 3100 blocks in a 
1.9 m by 1.9 m array. The depth of the Csl calorime- 
ter corresponds to 27 radiation lengths. Two 15 cm by 



15 cm holes are located near the center of the array for 
the passage of the two neutral beams. For electrons with 
energies between 2 and 60 GeV, the calorimeter energy 
resolution is below 1% and the nonlinearity is less than 
0.5%. The position resolution of the calorimeter is ap- 
proximately 1 mm. Downstream of the Csl calorimeter, 
there is a 10 cm lead wall, followed by a hodoscope used 
to reject hadrons hitting the calorimeter. 

The Kl — > 7r°e + e _ 7 decays were required to satisfy 
certain trigger requirements in order to be recorded. In 
particular, activity in a set of hodoseopes upstream of 
the Csl calorimeter had to be consistent with two tracks. 
Also, we required the event have at least one hit in one 
of the two upstream drift chambers. The event must 
deposit more than approximately 25 GeV of total en- 
ergy in the Csl calorimeter and no more than 0.5 GeV 
in the photon vetoes. The event is vetoed if it deposits 
more than 2.5 Mips in the hodoscope downstream of the 
calorimeter or more than 14 GeV in the vetos around the 
beam holes in the Csl calorimeter. The trigger includes 
a hardware cluster processor that counts the number of 
in-time calorimeter clusters of contiguous blocks of Csl 
with energies above 1 GeV[l5[. The total number of elec- 
tromagnetic clusters in the Csl calorimeter is required to 
be greater than or equal to four at the trigger level. 

After the events are read out, they must satisfy a soft- 
ware filter. This filter requires that each event have two 
charged tracks with a minimum of four clusters in the 
calorimeter. Each of the tracks must point to a clus- 
ter in the calorimeter and be consistent with an elec- 
tron hypothesis. The trigger requirements also select 
Kl — * ^°^d events where one neutral pion undergoes 
Dalitz decay, tt° — > e + e~7 (7r^,). These events are used 
for normalizing the Kl — > ir°e + e~j events, since their 
topology is very similar to that of our signal events. Be- 
cause of the similarity in topologies between the signal 
and normalization modes, many systematic effects can- 
cel. 



III. EVENT RECONSTRUCTION 

The offline analysis begins by requiring that each event 
have exactly two oppositely signed tracks and five in-time 
clusters with energies greater than 2.0 GeV, where an 
in-time cluster is one in which the cluster reconstructs 
to within 19 ns of the event time. The two tracks are 
required to point to two of the clusters and be consis- 
tent with a common decay vertex. From the three neu- 
tral clusters, we combine two to form the ir° candidate. 
There are three possible combinations and we choose the 
combination that reconstructs closest to the 7r° mass. 
Only events with a 77 invariant mass within 5 MeV/c 2 
of the nominal 7r° mass are accepted. The neutral de- 
cay distance vertex is used to determine the mass of the 
e + e~7 and e + e~777 combinations because it improves 
their mass resolution; the e + e~ tracks tend to be close 
together leading to poorer vertex resolution. The total 



3 



2%J> Kaon Mass 




0.40 0.45 0.50 0.55 0.60 

eera Mass GeV/c 2 



FIG. 2: The Kl — > it tt d invariant mass distribution for data 
(dots) and Kl — > 7r f V^ Monte Carlo (solid histogram). 



kaon energy, determined from the sum of cluster energies 
in the calorimeter, must lie between 30 and 210 GeV. 

To ensure that the two tracks are electrons, the recon- 
structed energy in the calorimeter divided by the mo- 
mentum determined by the spectrometer (E/p) of each 
track must be between 0.95 and 1.05. Backgrounds from 
K s decays and misreconstructed kaons can be reduced by 
requiring the decay vertex to reconstruct between 98 and 
157 meters downstream of the target, and the transverse 
momentum squared (p|n) for the event to be less than 
0.003 (GeV/c) 2 . The invariant mass for K L -> 7r°7r£, 
events is shown in Figure [2] The data and our Monte 
Carlo simulation agree quite well. 



IV. BACKGROUNDS TO K L -> ^e+e'-y 

After applying the above selection criteria, the remain- 
ing backgrounds consist mainly of Kl — > 7r°7r^ ) and 
K L — > 7r 7r°7r^, decays. The K L — > 7r°7r^, decays are 
more readily removed because the invariant masses of the 
e + e~7 and 77 combinations reconstruct around the mass 
of the 7T°. The majority of the Kl — > 7r°7r^, events are re- 
moved by requiring the reconstructed e + e~j mass of the 
best combination to be less than 0.110 GeV/c 2 or greater 
than 0.155 GeV/c 2 . However, when the wrong 77 com- 
bination is chosen, a restriction on the invariant masses 
is ineffective at reducing the background. Kl — > 7r°7r£, 
events can also contribute to the background if one of 
the final state particles is lost and is replaced by activ- 
ity in the detector that can mimic a final state particle. 
The Kl — > tt°tt°tt d events are more difficult to remove 
because we cannot use the same mass constraints as in 
the Kl — > 7r i/j case. However, kinematic and cluster 
shape variables have been developed to help to reduce 
the background to a manageable level. 

To remove misreconstructed Kl ^ir^ decays, we 
consider the two other possible 77 combinations. We 
take advantage of the correlations between the m 77 and 



m e + e - 7 distributions for these two combinations, form- 
ing a neural net from four variables. These four input 
variables are the reconstracted invariant 77 and e + e~7 
masses for each of the two remaining combinations. The 
neural net employed sixteen hidden nodes and was tuned 
on a sample of 27r^, and 7r°e + e~7 Monte Carlo. The out- 
put from the neural net ranges between zero and one. We 
rejected events where the neural net value was less than 
0.5. 

Backgrounds from Kl 7r 7r 7r^, come from two 
broad classes of events: events with missing photons and 
those with one or more photons that overlap or fuse to- 
gether in the Csl calorimeter. For events with missing 
photons, we use the photon vetoes to significantly reduce 
the amount of background. We require the maximum 
energy in any photon veto to be less than 0.1 GeV. To 
reduce backgrounds from events with overlapping pho- 
tons, we examined the calorimeter energies in a 3 x 3 ar- 
ray of crystals centered around the highest energy crystal 
of the cluster. For reference, a 3 x 3 array corresponds 
to approximately one Moliere radius. We compared these 
energies to energies from an ideal cluster shape and calcu- 
lated a photon shape \ 2 variable. This variable is shown 
in Figure [3] As can be seen, for the normalization mode, 
there is good agreement in this variable between the data 
and the Monte Carlo simulation. For the signal events, 
the background from Kl ~ > 7r°7r 7T£, events is signifi- 
cantly reduced by requiring a small value of shape % 2 - 
We require shape x 2 < 4. 

Kaon decays with missing photons will also exhibit a 
significant amount of missing energy when boosted to the 
center-of-mass. We take advantage of this effect by calcu- 
lating the longitudinal missing momentum in the center- 
of-mass (ppOkine). In the ppOkine versus m 777 plane, the 
signal events are well-separated from the Kl — > tt ^^ 
background. We define a two dimensional cut by employ- 
ing the following fourth-order polynomial: 

pp0kme maa , = A + B * (to 777 — xo) + C * (to 777 — xq) 2 
+ D * (to 777 — xq) 3 + E * (to 777 — Xo) 4 . 

where A = 3.9, B=-112.8, C=1256.6, D=-5861.8, 
E=10506.0 and x = 8.326 x 10~ 2 . The values of 
these parameters were chosen to maximize the signal-to- 
background ratio. Events with values of ppOkine greater 
than this value were rejected. This cut is shown in Fig. 2] 

Both Kl — > 7]" 7T£> and Kl — * 7r 7r°7i£, events can 
contribute to the final data sample if one of the elec- 
trons undergoes bremsstrahlung radiation. To reduce 
this background, we calculate the minimum distance be- 
tween the projection of the upstream segment of each 
electron to the Csl calorimeter and any photon clus- 
ter. Bremsstrahlung photons tend to have a small min- 
imum distance. As shown in Fig. [5J backgrounds from 
Kl — > 7r°7]/j and Kl — > 7r 7r 7i£, occur at small photon- 
to-track distances. To reduce this background, we require 
the minimum distance to be greater than 1.25 cm. 

Backgrounds due to external conversions of photons 
are negligible since we require the decay vertex to lie 
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FIG. 5: The minimum distance between the projection of 
the upstream electron track segment to the Csl calorimeter 
and any photon candidate. The dots represent the data and 
the solid histogram the sum of the background and signal 
Monte Carlo. The enhancement at low values of the minimum 
distance is the result of bremsstrahlung production. 



FIG. 3: The photon shape variable shape x 2 ■ The top plot 
shows the shape \ 2 variable for Kl — * t^°tvd with the dots rep- 
resenting the data and the solid histogram the Monte Carlo. 
In the bottom histogram the dots are the data after remov- 
ing the Kl — * tt ^ events, while the solid histogram shows 
the Kl — ► it°n°TY% Monte Carlo. The dashed histogram rep- 
resents the sum of the signal and background Monte Carlo 
samples. 



within the vacuum decay region, upstream of the vac- 
uum window. The requirement that the neutral and 
charged vertex be consistent with each other further re- 
duces the probability of external conversions contributing 
to the background since external conversions can only re- 
sult from interactions with material downstream of the 
vacuum window. 



V. RESULTS 
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FIG. 4: The ppOkine variable plotted versus the 777 invariant 
mass for a) Kl — ► ir -K% and b) Kl — ► 7r 7r 7r£> Monte Carlo 
events. The dark line represents the cut indicated in the text. 



After making these final selection criteria, we find the 
e + e~777 mass distributions shown in Figure [S] A clear 
peak at the kaon mass is seen, while the background is 
well-described by the sum of the 2it° D and 3^ back- 
ground Monte Carlo samples. We find a total of 139 can- 
didate events with an estimated background of 14.4 ± 2.5 
events. 

The K L — > 7r°e + e _ 7 branching fraction is determined 
from the following expression: 

B = (N 7r o e + e - 1 /N 27T o) X (e27r°/ e 7r°e+e-7) 
x B{K L -> 7r°7r°) x B(n° -> e+e'j) x 2. 

Here, N n o e + e - 1 represents the number of signal candi- 
dates, while N 2ll o represents the number of normalization 
events. The number of K L — > candidates is deter- 

mined by inverting the cut against Kl — > ^ir^ events 
and counting the number of events in the kaon mass re- 
gion from 0.490 to 0.510. In the above expresssion e 2 7r 
and e,rO e + e - 7 correspond to the reconstructed Kl — > 27r° 
and Kl — * 7r°e + e~7 acceptances, respectively. The 
factor of two occurs because there are two 7r° in each 
Kl — » 7r 7r£, event. In the previous analysis, the value of 
B(K L -> 7r°7r°) used was (9.36 ±0.2) x 10~ 4 . We are now 
using the most recent determination of B(Kl — > 7r°7r°) 
: -x.fi?) ± 0.08) x 10~ 4 . The value of B(tt° -> e+e'-f) used 
in both analyses is (1.198 ± 0.032) x 10" 2 . 
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FIG. 6: The e + e~777 invariant mass for events passing all se- 
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histogram represents the sum of the signal and background 
Monte Carlo. The background Monte Carlo is indicated by 
the shaded histogram. 



The acceptance for Kl — > "2^% events is 0.51% in the 
1997 data set and 0.61% in the 1999 data set. The differ- 
ence between the acceptances in the two data sets arises 
from the different magnetic fields used during the 1997 
and 1999 runs. We find 31,286 K L -> tt tt° d events in 
the 1997 data and 49,159 events in the 1999 data. This 
corresponds to a kaon flux of 2.88 x 10 11 and 3.78 x 10 11 
decays in the 1997 and 1999 data sets, respectively. The 
7r°e+e~7 acceptances are 0.90% and 1.02% for the 1997 
and 1999 data sets, respectively. These values are shown 
in Table [T] Using the numbers above, we obtain: 

B(A'i — ► 7r°e + e~7) = (1.49 ± 0.22) x 10" 8 (1997) 
B(K L -> 7r°e + e"7) = (1.77 ± 0.18) x 10" 8 (1999) 

The difference between the result in Ref. [l(| and our 
new 1997 measurement of the Kl — ► 7r°e + e _ 7 branch- 
ing ratio can be mainly attributed to the different values 
of ay used in determining the acceptance. The previous 
measurement used ay = —0.96 while our new analysis 
uses ay — —0.460. After accounting for this acceptance 
effect, we find that the two analyses are consistent with 
each other. Our new measurement supersedes the previ- 
ous KTeV measurement. 



Value 1997 1999 

Events in Data 47 92 

Background Events 2.7 11.7 

Normalization Events 31,286 49,159 
Signal Acceptance 0.91% 1.03% 

Normalization Acceptance 0.51% 0.61% 



TABLE I: Values used in branching ratio calculation. 



VI. SYSTEMATIC UNCERTAINTIES 

The largest systematic uncertainty results from the 
limited statistics in our background Monte Carlo sample. 
In total we generated approximately twice the statistics 
of the Kl ~ * 7r°7r 7i/j) data sample, and approximately 
three times the statistics of the normalization and sig- 
nal modes. This required generating about six billion 
Monte Carlo events. The next largest systematic uncer- 
tainty arises from the Kl and ir° — ► e + e~j branching 
ratios. The remaining effects can be broken down into 
two main classes: those that affect the background level 
and those that affect the signal or normalization accep- 
tance. The signal acceptance has a dependence upon the 
value of ay. This dependence can be characterized by 
B = 1.484 + 0.215 *a v + 0.1.170 * a\ + 0.417 * oi. We 
varied the value of ay between -0.41 and -0.51 Q, and 
found the acceptance changed by approximately 2.0%. 
The uncertainty in the background contributes approx- 
imately 0.5% to the total systematic uncertainty, and 
the remaining acceptance effects including the effects of 
apertures and cuts contribute about 0.6% to the total 
systematic error. All of the systematic errors are listed 
in Table HJ 



Systematic 


Error (%) 


MC Statistics 


4.2 


K L -> ttV and ir° D BR 


2.8 


av dependence 


2.0 


Signal acceptance 


0.6 


3ir ( £, and 2n% background 


0.5 


Total 


5.2 



TABLE II: Systematic uncertainties in percent. 

To obtain the final result, we took the weighted average 
of the 1997 and 1999 numbers, where we weighted by 
the statistical error. The systematic studies were done 
on the combined 1997 and 1999 analyses to take into 
account any correlations. Including the uncertainties due 
to the systematic effects, we find the following result: 
B(K L -► ^°e+e-7) = (1.62 ± 0.14 atat ± 0.09 syst ) x KT 8 . 

The relatively small value for the Kl 7r°e + e~7 
branching ratio means that this decay will not serve as a 
significant source of background to Kl — * 7r°e + e~. The 
7r°e + e~ mass distribution from Kl —* 7r°e + e~7 decays 
peaks around 0.250 GeV/c 2 , far enough away from the 
kaon mass that only a small fraction of the decays pose 
any risk of reconstructing in the Kl —> 7r°e + e~ signal 
region. 

To determine the value of ay from our data, we per- 
formed a maximum likelihood fit to the three Dalitz pa- 
rameters Z — m "n 9 , Q = ^fS and Y = — — ^ ££ . The 

variables E 1 and E ee are the energies of the photon and 
the e + e~ pair in the kaon CM frame, respectively. The 
value of ay that we obtained is —0.76 ± 0.16 ± 0.07. The 
major systematics are the background level and the sen- 
sitivity to the selection criteria. Our value for ay is con- 
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FIG. 7: The Dalitz variables used in the fit for ay: a) m e + e - 7 , 
b) Y and c) m e + e -. The dots are the data. The dashed his- 
togram is the sum of the background and signal Monte Carlo, 
and the filled histogram is the background Monte Carlo. 



sistent with the recent published values but our errors 
are significantly larger. The distributions of the Dalitz 
variables are shown in Fig. \7\ We find good agreement 
between the data and the Monte Carlo generated with 
the central ay value. The shapes of the m e + e - 7 and Y 
variables are similar to those seen in Kl — > 7r°77 de- 
cay. In principle the e + e _ mass distribution should rise 
sharply near threshold. However, due to the detector 
acceptance this peaking is suppressed. 

In Ref.[l| the 0(p 6 ) calculation predicts the Kl — > 
7r°e + e - 7 branching ratio to be 2.4 x 10~ 8 . This pre- 
diction depends upon the value of ay. For a value of 
ay = —0.46 the branching ratio prediction is 1.51 x 10~ 8 , 
which is consistent with our measurement. The 0(p i ) 
calculation, meanwhile, predicts a branching ratio of 
1 x 10 -8 . Therefore, our latest measurement disfavors 



tion. 



VII. CONCLUSIONS 



We have determined the branching ratio B(Kl —> 
7r°e + e~7) using the combined 1997 and 1999 data sets 
from the KTeV experiment. The statistics represents 
a factor of 2.5 over our published 1997 result. Com- 
pared to our previous result, this analysis utilizes a 
number of new analysis techniques and employs an im- 
proved understanding of the backgrounds. We deter- 
mine the branching ratio to be ~Q(Kt ~ > 7r°e + e _ 7) = 
(1.62±0.14 stat ±0.09 syst ) x 1CT 8 . The most recent mea- 
surement of ay [|| suggests that the decay Kl — > 7r°e + e~ 
is dominated by a CP violating amplitude. While the 
statistics are low, the value of ay from our determination 
is consistent with these conclusions. Our branching ratio 
measurement confirms the value of ay measured by NA48 
and indicates that 0(p 6 ) terms are important for mod- 
eling this decay mode. The measured branching ratio is 
within 0.7 a of the 0(p 6 ) prediction and about 3.7 a from 
the 0(p 4 ) prediction. A factor of twenty increase in the 
statistics of Kl — > 7r°e + e~7 would make the measure- 
ment of ay competitive with the current best measure- 
ments. Because of the kinematics and small branching 
ratio for Kl — > 7r°e + e _ 7, this decay will not constitute 
a large background in future searches for Kl — > 7r°e + e~. 
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